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The study of exclusive processes in the future electron-ion (eA) colliders will be an important tool 
to investigate the QCD dynamics at high energies as they are in general driven by the gluon content 
of the target which is strongly subject to parton saturation effects. In this paper we compute the 
coherent and incoherent cross sections for the deeply virtual Compton scattering (DVCS) process 
relying on the color dipole approach and considering different models for the dipole - proton scat¬ 
tering amplitude. The dependencies of the cross sections with the energy, photon virtuality, nuclear 
mass number and squared momentum transfer are analysed in detail. We demonstrate that the ratio 
between the incoherent and coherent cross sections decreases at smaller values of Q 2 and increases 
at smaller values of A. We show that the coherent cross section dominates at small t and exhibits 
the typical diffractive pattern, with the number of dips in the range |t| < 0.3 GeV 2 increasing with 
the mass atomic number. Our results indicate that the position of the dips are independent of the 
model used to treat the dipole - proton interaction as well as of the center-of-mass energy. 

PACS numbers: 12.38.-t,13.60.Hb, 24.85.+P 


I. INTRODUCTION 

The study of the hadronic structure in the non-linear 
regime of the Quantum Chromodynamics (QCD) is one 
of the main goals of the future Electron - Ion Collid¬ 
ers (EIC) [ll-(|. It is expected that the resulting exper¬ 
imental data will be able to determine the presence of 
gluon saturation effects, the magnitude of the associated 
non-linear corrections and what is the correct theoreti¬ 
cal framework for their description. This expectation is 
easily understood, if we analyse the Bjorken-a; and nu¬ 
clear mass number A behaviours of the nuclear satura¬ 
tion scale, Q s ,a , which determines the onset of non-linear 
effects in the QCD dynamics Q. Assuming that the satu¬ 
ration scale of the nucleus is enhanced with respect to the 
nucleon one by the oomph factor A 3 , it can be expressed 
by the parametrization Q 2 A = As xQl (yf ) A - As a con¬ 
sequence nuclei are an efficient amplifier of the non-linear 
effects. In Fig. [Tjwe present the theoretical expectations 
for the saturation scale as a function of x and A , consid¬ 
ering the parameters Q\ = 1.0 GeV 2 , xq = 1.632 x 10 -5 
and A = 0.2197 as in Ref. We can observe that, 
while in the proton case we need very small values of x 
to obtain large values of Q 2 , in the nuclear case a similar 
value can be obtained for values of x approximately two 
orders of magnitude greater. Consequently, the parton 
density that is accessed in electron - nucleus (eA) colli¬ 
sions would be equivalent to that obtained in an electron 
- proton collider at energies that are at least one order of 
magnitude higher than at HERA. 

The enhancement of the non-linear effects with the nu¬ 
clear mass number has motivated the development of an 
intense phenomenology about the implications of these 
effects in inclusive and diffractive observables which could 
be measured in the Electron - Ion Colliders [3-ill. These 
studies indicate that the analysis of inclusive observables, 


as e.g. the nuclear structure functions and F£, prob¬ 
ably will not be the best way to obtain a signature of the 
non-linear effects due to the large uncertainty present in 
the collincar predictions at small x. In contrast, diffrac¬ 
tive events in eA collisions are expected to be a smoking 
gun for the gluon saturation physics. Basically, satura¬ 
tion physics predicts that in the asymptotic limit of very 
high energies, diffractive events will be half of the total 
cross section, the other half being formed by all inelas¬ 
tic processes UM- In the kinematical range probed by 
the future EIC, diffractive processes are expected to con¬ 
tribute by a large amount (~ 20%) to the total cross sec¬ 
tion 0. This observation have motivated more detailed 
studies of the diffractive interactions, with special atten¬ 
tion to exclusive diffractive vector meson production and 
the deeply virtual Compton scattering (DVCS), which 
are experimentally clean and can be unambiguously iden¬ 
tified by the presence of a rapidity gap in final state. As 
these exclusive processes are driven by the gluon con¬ 
tent of the target, with the cross sections being propor¬ 
tional to the square of the scattering amplitude, they are 
strongly sensitive to the underlying QCD dynamics. In 
particular, in Refs. [u| HU we have presented a system¬ 
atic analysis of these processes in terms of the non-linear 
QCD dynamics and derived the energy dependence of 
the total vector meson {V = p, <j> and J/iH) and DVCS 
cross sections. A shortcoming of these studies is that the 
distributions on the square of the momentum transfer 
(t), which are an important source of information about 
the spatial distribution of the gluons in a nucleus and 
about fluctuations of the nuclear color fields, were not 
estimated. Recently, these distributions were calculated 
in Refs. 0,ii, where a comprehensive study about 
the vector meson production in eA colliders considering 
different assumptions were presented. Our goal in this 
paper is to complement these studies by presenting a de- 
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FIG. 1: (Color online) Nuclear saturation scale as a function 
of the Bjorken x and nuclear mass number A. 


tailed analysis of the nuclear DVCS (For previous studies 
see, e.g. Refs. [23l - l28j ). In particular, we update the pre¬ 
vious calculations considering the more recent models for 
the dipole - proton scattering amplitude and presenting 
the predictions for the f-dependence of the nuclear DVCS 
cross section. 

This paper is organized as follows. In the next Sec¬ 
tion we present a brief review of the description of the 
nuclear DVCS in the color dipole formalism. In Section 
m we discuss the models for the dipole - proton scat¬ 
tering amplitude used as input in our calculations and 
present a comparison between its predictions. In Sec¬ 
tion m we present our results for the energy, virtuality 
and momentum transfer dependencies of the DVCS cross 
sections. Finally, in Section [V] we summarize our main 
conclusions. 


where 

(A f A ) = J d 2 r j dz%(r,z)Rf A {x,r,b) V^{r,z,Q 2 )( 2 ) 


and J\f A (x,r,b) is the forward dipole-nucleus scattering 
amplitude for a dipole with size r and impact parameter 
b which encodes all the information about the hadronic 
scattering, and thus about the non-linear and quantum 
effects in the hadron wave function. As in our previous 
studies [H, HU , in what follows we will use in our calcu¬ 
lations the model proposed in Ref. [30], which describes 
the current experimental data on the nuclear structure 
function as well as includes the impact parameter de¬ 
pendence in the dipole nucleus cross section (For details 
see Ref. 0 )- In this model the forward dipole-nucleus 
amplitude is given by 


M A {x, r, b) = 1 — exp 



a dp (x,r 2 )AT A (b) 


,(3) 


where a dp is the dipole-proton cross section and T A {b) 
is the nuclear profile function, which is obtained from 
a 3-parameter Fermi distribution for the nuclear density 
normalized to 1. This equation, based on the Glauber- 
Gribov formalism 0, sums up all the multiple elastic 
rescattering diagrams of the qq pair and is justified for 
large coherence length, where the transverse separation 
r of partons in the multiparton Fock state of the pho¬ 
ton becomes a conserved quantity, i.e. the size of the 
pair r becomes eigenvalue of the scattering matrix. The 
differential distribution with respect to the squared mo¬ 
mentum transfer t for coherent interactions is given by 


dtJcoh 

dt 

where 


1 

167T 


A^ A ^ A (x,Q 2 ,t) 


(4) 


A 1 a_> 7 j 4 (x, Q 2 , t) = 2 J d 2 be l ^' A {U A {x, r, 6 )),(5) 


II. NUCLEAR DVCS IN THE COLOR DIPOLE 
PICTURE 

Let us start presenting a brief review of the deeply 
virtual Compton scattering (DVCS) in electron - ion col¬ 
lisions. In the color dipole approach the exclusive pro¬ 
duction 7 * A —»• 'yY in electron-nucleus interactions at 
high energies can be factorized in terms of the fluctu¬ 
ation of the virtual photon into a qq color dipole, the 
dipole-nucleus scattering by a color singlet exchange and 
the recombination into the exclusive final state 7 . This 
process is characterized by a rapidity gap in the final 
state. If the nucleus scatters elastically, Y = A, the pro¬ 
cess is called coherent production and can be represented 
by the diagram in Fig. [2] (upper panel). The correspond¬ 
ing integrated cross section is given in the h igh energy 
regime (large coherence length: l c R A ) by [13,0 

a coh (7*A ->■ 1 A) = J d 2 b (U A (x, r, b)) 2 (1) 


with A 2 = —t and N A given by Eq. ([3]). 

On the other hand, if the nucleus scatters inclastically, 
i.e. breaks up (Y = A'), the process is denoted inco¬ 
herent production and can be represented as in Fig. [2] 
(lower panel). In this case one sums over all final states 
of the target nucleus, except those that contain particle 
production. Therefore we have: 


<r inc (7*A -a 7 A') = 


| ImA(s, t = 0 )| 2 
16t tB 


( 6 ) 


where at high energies (l c R A ) p 
\lmA\ 2 = J d 2 b AT A (b) l^a dp exp 


-- <J dp AT A (b ) 


(7) 


with the t slope B being the same as in the case of a 
nucleon target. In the incoherent case, the qq pair atten¬ 
uates with a constant absorption cross section, as in the 
Glauber model, except that the whole exponential is av¬ 
eraged rather than just the cross section in the exponent. 
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For the calculation of the differential cross section da/dt 
for incoherent interactions we apply for the DVCS pro¬ 
cess the treatment presented in Ref. [r|. Consequently, 
we have that 


da i nc 
dt 


±-((\A^ A ~^ A '(x,Q 2 ,t)\ 2 )), 


( 8 ) 


FIG. 2: Deeply Virtual Compton Scattering in coherent (up¬ 
per panel) and incoherent (lower panel) eA collisions. 


with the scattering amplitude A 7 A ^ A being approxi¬ 
mated by fl9j ] 


J 


|^ 7 * A ^ 7 A '( X ; q 2^)|2 = 167 r 2 B 2 J fpte - B , A a j ^( Xir ) J ^( x y ) ATA ( b ) 

x exp | — 27r(A — l)BpTA(b) [J\f p (x, r) + J\f p (x, r 7 )] (9) 


where B p is associated to the impact parameter profile 
function in the proton and J\T p (x , r) is the dipole - proton 
scattering amplitude to be discussed in more detail in the 
next Section. 

The functions \R 7 . ( 2 :, r) and ^tj{z, r) in the Eq. © 
are the light-cone wave functions of the virtual photon 
in the initial state and the real photon in the final state, 
respectively. The variable r defines the relative trans¬ 
verse separation of the pair (dipole) and z (1 — z) is 
the longitudinal momentum fraction of the quark (an¬ 


tiquark). In the dipole formalism, the light-cone wave 
functions ^(z, r) in the mixed representation (r, z) are 
obtained through a two dimensional Fourier transform of 
the momentum space light-cone wave functions \l/( z, k). 
The photon wave functions are well known in literature 
!3'J . In the DVCS case, as one has a real photon at the 
final state, only the transversely polarized overlap func¬ 
tion contributes to the cross section. Summed over the 
quark helicities, for a given quark flavour / it is given by 

m, 


j 




Nc Q!em6 f 
2 ^ 


{ [z 2 + x 2 ] £iK 1 (£ 1 r)£ 2 Ki(s 2 r) + m 2 K 0 (eir)K 0 (£ 2 r)} , 


( 10 ) 


where we have defined the quantities e 2 2 = zz Q\ 2 + m 2 
and 5 = (1 — z). Accordingly, the photon virtualities 
are Q\ = Q 2 (incoming virtual photon) and Q\ = 0 
(outgoing real photon). 


III. QCD DYNAMICS 

The DVCS cross sections in eA collisions are expressed 
in terms of the dipole-proton cross section ad pi which in 
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the eikonal approximation is given by: 

a dp (x,r) = 2 J d 2 bN p (x,r,b) , ( 11 ) 

where N p (x,r , 6 ) is the imaginary part of the forward 
amplitude for the scattering between a small dipole (a 
colorless quark-antiquark pair) and a dense hadron tar¬ 
get, at a given rapidity interval Y = ln(l/a;). The dipole 
has transverse size given by the vector r = x — y, where 
x and y are the transverse vectors for the quark and anti¬ 
quark, respectively, and impact parameter b = (x + y)/2. 
At high energies the evolution with the rapidity Y of 
M p (x , r, ft) is g iven in the Color Glass Condensate (CGC) 
formalism [331 ] by the infinite hierarchy of equations, the 
so called Balitsky-JIMWLK equations [33|,(34j], which re¬ 
duces in the mean field approximation to the Balitsky- 
Kovchegov (BK) equation [H, [HJ (For a detailed dis¬ 
cussion about the subject see, e.g., Ref. [3fij|). In re¬ 
cent years, the running coupling corrections to BK evo¬ 
lution kernel was explicitly calculated [§ 3 , S], includ¬ 
ing the a s Nf corrections to the kernel to all orders, 
and its solution studied in detail [39j, [4(J. Basically, 
one has that the running of the coupling reduces the 
speed of the evolution to values compatible with exper¬ 
imental ep HERA data [4lU43l|. The numerical solu¬ 
tions of the running coupling BK equation presented in 
Refs. [4ll - [43[ assumed the translational invariance ap¬ 
proximation, which implies J\f p (x,r,b) = J\T p (x,r)S(b) 
and ad,p(x,r ) = <to • Af p (x,r ), with the normalization of 
the dipole cross section (er 0 ) being fitted to data. It is 
important to emphasize that the normalization uq de¬ 
fines Bp in Eq. ©, since they are related: B p = cj 0 /At:. 
Moreover, J\f p (x,r ) is the main input for the calculation 
of the incoherent DVCS cross section [See Eq. ©]. In 
what follows we will consider the solutions for the scat¬ 
tering amplitude obtained in Refs. [4l|, |43|, denoted re¬ 
spectively by rcBK and AAMQS, which differ in the set 
of experimental data used to constrain the initial condi¬ 
tions and in the treatment of the heavy quark contribu¬ 
tions. While in the rcBK model [4l| the experimental 
data for the inclusive and longitudinal proton structure 
functions measured in ep collisions at HERA were con¬ 
sidered in the fit, the more precise HERA data for the re¬ 
duced cross sections were used in Ref. (H}. Moreover, in 
the rcBK model the contribution of the heavy quarks was 


disregarded. In contrast, in Ref. [43[ the authors have 
performed a non-linear QCD analysis of new HERA data 
at small-a; including heavy quarks. We denote the corre¬ 
sponding solution of the running coupling BK equation 
by AAMQS(H). For comparison, we also will consider the 
solution derived in Ref. {43j, denoted AAMQS(L) here¬ 
after, which was obtained taking into account only light 
quarks (See Ref. [43| for more details). 

A shortcoming of the solutions presented in Refs. [H- 
fjdj is that they disregard the impact parameter depen¬ 
dence. Unfortunately, impact-parameter dependent nu¬ 
merical solutions to the BK equation are very difficult 
to obtain (44|. Moreover, the choice of the impact- 
parameter profile of the dipole amplitude entails intrinsi¬ 
cally nonperturbative physics, which is beyond the QCD 
weak coupling approach of the BK equation. In fact, the 
BK equation generates a power law Coulomb-like tail, 
which is not confining at large distances and therefore 
can violate the Froissart bound [45[ . As demonstrated in 
Refs. [44] . [ifij , confinement can be modelled by the mod¬ 
ification of the dipole evolution kernel of the BK equa¬ 
tion by including an effective gluon mass, which implies 
that the Froissart bound limit is satisfied. An alterna¬ 
tive is to construct phenomenological models which in¬ 
corporate the expected 6 -dependence of the scattering 
amplitude. Although a complete analytical solution of 
the BK equation is still lacking, its main properties at 
fixed b are known: (a) for the interaction of a small 
dipole (r = |r| <C 1 /Q s ), Af p (x,r,b) ss r 2 , implying 
that this system is weakly interacting; (b) for a large 
dipole (r 1 /Q s ), the system is strongly absorbed and 
therefore Af p (x,r,b) ps 1. This property is associated 
to the large density of saturated gluons in the hadron 
wave function. In the last years, several groups have 
constructed phenomenological models which satisfy the 
asymptotic behaviours of the BK equation in order to fit 
the HERA and RHIC data (See e.g. Refs. [13, S3-E3l)- 
In particular, in Ref. [32| the authors have proposed a 
generalization of Iancu - Itakura - Munier model jUj], 
which is based on the CGC formalism, with the inclusion 
of the impact parameter dependence in the dipole - pro¬ 
ton scattering amplitude (denoted bCGC hereafter). In 
the bCGC model the dipole - proton scattering amplitude 
is given by |42 


A f p (x, r , b) 


A/o 

1 — exp -A 


^ 2 ( 7a + ln(2 ( r ^ (l>)) 

In 2 ( BrQ B (b )) 


) 


rQ s (b) < 2 
rQ s {b) > 2 


( 12 ) 


with k = x" (is)/ x' (is ), where x is the LO BFKL char¬ 
acteristic function. The coefficients A and B are deter¬ 
mined uniquely from the condition that M p (x , r*, fo), and 


its derivative with respect to rQ s (b), are continuous at 
rQ s (b ) = 2. In this model, the proton saturation scale 
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Q s (b) depends on the impact parameter: 


Qs O) = Qs{x,b) = 



exp 


AVll 

2-Bcgc /_ 


(13) 


The parameter Bcgc was adjusted to give a good de¬ 
scription of the f-dependence of exclusive J/ij) photopro¬ 
duction. Moreover, the factors A/q and % were taken 
to be free. In this way a very good description of F 2 
data was obtained. Recently, this model has been im¬ 
proved by the fitting of its free parameters considering 
the high precision combined HERA data j53|. The set 
of parameters which will be used here are the following: 

= 0.6599, k = 9.9, B CG c = 5.5 GeV" 2 , Af 0 = 0.3358, 
Xq = 0.00105 and A = 0.2063. In order to estimate 
the model dependence of our predictions, in what fol¬ 
lows we also will present the predictions obtained con¬ 
sidering the phenomenological models discussed in Refs. 
0 S 3 ; which assume the factorization of the impact 
parameter dependence and propose different forms for 
A f p (x, r). In particular, the GBW model (13 assume that 
Af p (x,r) = 1 — and the IIMS model 0, HU 

can be obtained from Eq. m disregarding the impact 
parameter dependence of the saturation scale. For the 
GBW model we use the following set of parameters HU : 
Ib = 3x 10~ 4 , A = 0.288 and cjq = 23.03 mb. For the 
IIMS model we assume that (U A/q = 0.7, = 0.7376, 

k = 9.9, A = 0.2197, x 0 = 1.632 x 10“ 5 and a 0 = 27.28 
mb. 

Before presenting our results for the nuclear DVCS 
cross sections in the next Section, we compare the 
dipole - proton scattering amplitudes predicted by the 
solutions of the running coupling BK equation (rcBK 
and AAMQS) with those from phenomenological models 
(IIMS and GBW). In Fig. [3] we analyse the pair separa¬ 
tion dependence of the dipole scattering amplitude A f p 
for two different values of x. For comparison we also 
present the predictions of the bCGC for two distinct val¬ 
ues of the impact parameter. It is important to empha¬ 
size that in Ref. [43j , when the heavy quark contribution 
is included in the analysis, the authors have considered 
different initial values of the parameters in the initial con¬ 
dition for light and heavy quarks. It implies two different 
parametrizations for the AAMQS (H) solution, one asso¬ 
ciated to light and other to heavy quarks. It explains 
the presence of two AAMQS(H) curves in Fig. [3] We 
observe that while the GBW and IIM parametrizations 
present a similar behaviour for small r 2 , the rcBK and 
AAMQS one predict a smoother dependence. These dif¬ 
ferences are associated to the distinct behaviours of the 
linear regime assumed by the models. In the saturation 
regime (large r 2 ), the GBW and IIM parametrizations 
saturates for large pair separations, while the rcBK and 
AAMQS one still present a residual dependence, demon¬ 
strating that the asymptotic regime is only reached for 
very large pair separations. The delayed saturation of 
the scattering amplitude at small values of x is still more 
pronounced in the bCGC predictions. The characteristic 



r 2 [GeV' 2 ] 


FIG. 3: (Color online) Dependence of the dipole - proton scat¬ 
tering amplitude in the squared pair separation r 2 at different 
values of x. 


feature which is evident in the GBW and IIMS models is 
that the dipole scattering amplitude saturates for smaller 
dipoles when x assumes smaller values, which is directly 
associated to the energy dependence of the proton sat¬ 
uration scale, which defines the onset of the saturation 
regime. Finally, the different models presented in Fig. [3] 
predict very distinct behaviours for the transition region 
between small and large pair separations. As we will 
demonstrate in the next Section, these differences has 
strong influence on the energy behaviour of the DVCS 
cross section. 


IV. RESULTS 


We initially estimate the energy dependence of the 
DVCS cross section in ep collisions. As in Ref. [50 we 
take into account the corrections associated to the skew¬ 
ness factor R g and to the real part of the scattering ampli¬ 
tude. Basically, the expressions for the DVCS cross sec¬ 
tions presented in Section [TT] have been multiplied by the 
factor Rg(l+ /3 2 ), with the factor R g taking into account 
that the gluons emitted from the quark and antiquark 
into the dipole carry different momentum fractions, and 
the factor /3 being the ratio of real to imaginary parts 
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FIG. 4: (Color online) Energy dependence of the DVCS cross 
section in ep collisions for two values of the photon virtuality 
Q 2 . Data from HI Collaboration J54j. 


of the scattering amplitude (For details see, e.g. Ref. 
0 )- Moreover, when using as input in our ep calcula¬ 
tions the dipole models which disregard the impact pa¬ 
rameter dependence (rcBK, AAMQS, GBW and IIMS), 
we will assume that & = -g • ^-|t=o> with the slope param¬ 
eter B being given by the experimental parametrization 
proposed in Ref. [54]: B (Q 2 ) = a[l — 61og(Q 2 / Qq)\, 
with a = 6.98 ± 0.54 GeV 2 , b = 0.12 ± 0.03 and Qg = 2 
GeV 2 . In contrast, in the case of the bCGC model, the 
cross section will be calculated by the integration of the 
differential cross section 4? as given in Eq. (4) of the 
Ref. jH]]. In Fig. [I] we present our predictions for the 
DVCS cross section in ep collisions as a function of the 
photon-target c.m.s energy, W, for two values of the pho¬ 
ton virtuality: Q 2 = 8 and 25 GeV 2 . For comparison 
the HERA data from the HI Collaboration [HiJ are also 
presented. We have that the GBW predictions are an 
upper bound for the cross section, while the bCGC pre¬ 
dictions can be considered a lower bound. The rcBK and 
AAMQS(L) predictions are similar, which is expected, 
since they are based on similar assumptions. In contrast, 
the AAMQS (H) prediction is similar to the IIMS and 
bCGC predictions. In the HERA kinematical regime, 
the difference between the predictions is not large, but 
becomes a factor 2 in the range which will be probed by 



FIG. 5: (Color online) Energy dependence of the coherent 
(upper panel) and incoherent (lower panel) DVCS cross sec¬ 
tions for A = 2 (D) and Q 2 = 8 GeV 2 . 


the future ep colliders Q. This huge difference is directly 
associated to the distinct behaviours for the dipole scat¬ 
tering amplitude in the transition region between small 
and large dipoles observed in Fig. [3] We have that the de¬ 
scription of the experimental data by the different models 
is satisfactory, in particular at Q 2 = 8 GeV 2 , where we 
expect a larger contribution of the saturation effects. In 
contrast, for large values of Q 2 , the main contribution 
comes from smaller dipoles, which implies that we are 
probing the linear regime. In this regime, effects associ¬ 
ated to the evolution in Q 2 , which are not included in the 
dipole models considered in this paper, can contribute. 
Consequently, in what follows we will restrict our analy¬ 
sis of DVCS cross sections to values of Q 2 smaller than 
15 GeV 2 . 

Lets now estimate the coherent and incoherent DVCS 
cross sections for different values of the nuclear mass 
number: A = 2 (D), A = 40 (Ca) and 208 (Pb). These 
results are an update of those presented in Ref. EH cor¬ 
responding to AAMQS predictions and light nuclei. For 
A = 2, the oomph factor As is close to one, implying a 
small enhancement of non-linear effects in comparison to 
ep collisions. However, the study of coherent and incoher¬ 
ent interactions in eD collisions can be useful to test the 
color dipole formalism, as well as to compare its predic- 
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W [GeV] 

FIG. 6: (Color online) Energy dependence of the coherent 
DVCS cross section for two different nuclei and Q 2 = 8 GeV 2 . 


tions with those obtained using very distinct frameworks 
(See, e.g. Refs. For completeness, in Fig. [5] we 

present the predictions for the energy dependence of the 
coherent and incoherent eD cross sections considering as 
input the different models for the dipole - proton scat¬ 
tering amplitude discussed before. On the other hand, 
for A = 40 (Ca) and 208 (Pb) we expect the enhance¬ 
ment of the non-linear effects, which makes the study of 
the DVCS process useful to constrain the magnitude of 
these effects. Consequently, considering the main goal of 
this paper, in what follows we will restrict our analysis 
to these two nuclei. In Fig. [6] we present our predictions 
for the energy dependence of the coherent cross sections 
for a fixed virtuality and in Fig. [7] the Q 2 dependence 
for a fixed energy, W = 100 GeV. Our results for A = 40 
and 208 agree with the bCGC and rcBK predictions pre¬ 
sented in Ref. [18]. We obtain that the cross sections 
increase with the energy and nuclear mass number and 
decrease with the photon virtuality. As expected from 
Fig. m the rcBK and AAMQS(L) predictions are very 
similar in the energy range analysed (W < 200 GeV). 
At Q 2 = 8 GeV the bCGC, IIMS and AAMQS(H) are 
almost identical. However, as can be verified in the Fig. 
0 its predictions differ at smaller values of Q 2 , due to 
the different behaviour predict by these models in the 
transition region between small and large dipoles. In the 



FIG. 7: (Color online) Dependence on the photon virtuality 
Q 2 of the coherent DVCS cross section for two different nuclei 
and W = 100 GeV. 


energy range which will be probed in future eA colliders 
([W < 140 GeV), the two AAMQS predictions differ by 
ss 50%. These conclusions are also valid for the incoher¬ 
ent cross sections presented in Figs. [8] and [9] In Fig. [TO] 
we present the energy dependence of the ratio between 
the incoherent and coherent cross sections for different 
nuclei and Q 2 using the rcBK model. We obtain that the 
ratio decreases at smaller values of Q 2 and increases at 
smaller values of A. It is important to emphasize that we 
have verified that our predictions are almost independent 
of the model used for the dipole - proton scattering am¬ 
plitude. Finally, we have obtained that the A dependence 
of our predictions can be approximated by a co h oc A 1 ' 3 
and (Ji nc oc A 0 75 for Q 2 = 8 GeV 2 and W = 120 GeV. 

Lets now extend our analysis for the differential distri¬ 
butions da/dt for coherent and incoherent interactions. 
In Fig. [FT] we present our predictions for the coherent 
DVCS cross section for different nuclei and Q 2 = 8 GeV 2 . 
We obtain that the coherent cross section clearly exhibits 
the typical diffractive pattern, with the dips in the range 
|i| < 0.3 GeV 2 increasing with the mass atomic num¬ 
ber. Moreover, the positions of the dips are almost in¬ 
dependent of the dipole - proton model used as input in 
the calculations. In Fig. [T2] we present our predictions 
for the incoherent cross section. In this case no dip is 



















































W [GeV] 

FIG. 8: (Color online) Energy dependence of the incoherent 
DVCS cross section for two different nuclei and Q 2 = 8 GeV 2 . 


present and the predictions of the different models being 
similar at large t. A more detailed comparison between 
the coherent and incoherent cross sections is presented 
in Fig. [13l where we consider distinct values of the en¬ 
ergy (upper panel) and different nuclei (lower panel). As 
expected from the studies of the exclusive vector meson 
production [Sill Hi , we obtain that coherent DVCS 
cross section dominates at small t (|f| • R 2 A /3 Cl), the 
signature being a sharp forward diffraction peak. On the 
other hand, incoherent DVCS production dominates at 
large t (\t\-R\/3 1), the t-dependence being to a good 

accuracy the same as in the production off free nucleons. 
We obtain that the magnitude of the differential cross 
sections at |f| = 0 increases with the energy and the nu¬ 
clear mass number, as expected from the analysis of the 
total cross sections, and that the positions of the dips are 
almost energy independent. 


V. CONCLUSIONS 

The study of the hadronic structure using hard exclu¬ 
sive reactions in deep inelastic scattering (DIS) processes, 
such as deeply virtual Compton scattering (DVCS) and 
vector meson production, have been one of the main fo¬ 
cuses of hadronic physics in the last years. These pro¬ 



FIG. 9: (Color online) Dependence on the photon virtuality 
Q 2 of the coherent DVCS cross section for two different nuclei 
and W = 100 GeV. 


cesses have been the subject of intensive theoretical and 
experimental investigations. One of the main motivations 
for these studies is the possibility to probe the QCD dy¬ 
namics at high energies, driven by the gluon content of 
the target (proton or nucleus) which is strongly subject 
to non-linear effects (parton saturation) effects. As in eA 
collisions the coherent contributions from many nucleons 
effectively amplify the gluon density probed, the study of 
exclusive observables in the future electron - ion colliders 
is the ideal scenario to constrain the QCD dynamics at 
high energies. In this paper we have presented a detailed 
study of the nuclear DVCS process, which complement 
the existing analysis about vector meson production in 
eA collisions. Although the DVCS cross section is smaller 
than the vector meson one, it is not affected by the the¬ 
oretical uncertainties associated to the scarce knowledge 
of the vector meson wave functions. Consequently, the 
DVCS process can be considered a direct probe of the 
QCD dynamics which describes the dipole - target inter¬ 
action. We have assumed in our studies that the dipole 
- nucleus scattering amplitude can be described in terms 
of the dipole - proton one and considered different mod¬ 
els for J\f p (x,r). In particular, we have extended the 
results obtained in Ref. [jj|, presenting the first time 
the predictions of the AAMQS model for the DVCS pro- 
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tualities. We have verified that the incoherent contribu¬ 
tion increases with the photon virtuality and decreases 
with the nuclear mass number. Moreover, we estimated 
the differential cross sections for the nuclear DVCS and 
have analysed in detail the dependence on the squared 
momentum transfer. As for vector meson production, 
we have obtained that the coherent cross section domi¬ 
nates at small t and the incoherent one at large t. We 
have demonstrated that the number of dips at small t 
increases with the atomic number, with the position of 
the dips being almost independent of the model used to 
treat the dipole - proton interaction and the center-of- 
mass energy. Both results are robust predictions from 
the saturation physics, which can be used to challenge 
the treatment of the non-linear QCD dynamics in the 
kinematical range of future electron - ion experiments. 
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